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ABSTRACT 

There is much debate on how high- mass star formation varies with environment, and whether the 
sparsest star-forming environments are capable of forming massive stars. To address this issue, we 
have observed eight apparently isolated OB stars in the SMC using HSTs Advanced Camera for 
Surveys. Five of these objects appear as isolated stars, two of which are confirmed to be runaways. 
The remaining three objects are found to exist in sparse clusters, with < 10 companion stars revealed, 
having masses of 1 - 4 M . Stochastic effects dominate in these sparse clusters, so we perform Monte 
Carlo simulations to explore how our observations fit within the framework of empirical, galactic 
cluster properties. We generate clusters using a simplistic -2 power-law distribution for either the 
number of stars per cluster (AT*) or cluster mass (M c \). These clusters are then populated with 
stars randomly chosen from a Kroupa IMF. We find that simulations with cluster lower- mass limits of 
M c \,\o > 20-M© and A^io > 40 match best with observations of SMC and Galactic OB star populations. 
We examine the mass ratio of the second-most massive and most massive stars m max ^/romax, finding 
that our observations all exist below the 20th percentile of our simulated clusters. However, all of our 
observed clusters lie within the parameter space spanned by the simulated clusters, although some 
are in the lowest 5th percentile frequency. These results suggest that clusters are built stochastically 
by randomly sampling stars from a universal IMF with a fixed stellar upper-mass limit. In particular, 
we see no evidence to suggest a rn ma x — M c \ relation. Our results may be more consistent with core 
accretion models of star formation than with competitive accretion models, and they are inconsistent 
with the proposed steepening of the integrated galaxy IMF (IGIMF) . 

Subject headings: galaxies: Magellanic Clouds - galaxies: star clusters - galaxies: stellar content - 
stars: early-type - stars: formation 



1. INTRODUCTION 

t> ' Most observational properties of galaxies and stellar 
C*~) , populations are directly influenced by star formation. 
C ■ Massive stars, although small in number, disproportion- 
ally affect observables such as the integrated light from 
*n | galaxies, feedback effects, star formation rates, and many 
f^i • others. However, there is significant debate regarding the 
— ' conditions under which massive stars form. The compet- 
\ itive accretion model of star formation requires that a 
7-H . population of low-mass stars must form in the presence 
!LJ 1 of a high- mass star (e.g. Bonnell et al. 2004), while the 
. £^ | core accretion model of star formation allows for massive 
stars to form in relative isolation (e.g. Krumholz et al. 
2009). Empirical studies are similarly divided on the coil- 
ed ' nection between cluster mass and massive star formation. 
A physical relationship between the two (e.g. Weidner & 
Kroupa 2006, hereafter WK06) would indicate that mas- 
sive stars always form in clusters, while random statisti- 
cal sampling (e.g. Elmegreen 2000) would indicate that 
massive stars preferentially form in clusters. One obser- 
vational method to differentiate between these theories 
is to examine the sparsest environments where massive 
stars are found. A targeted study of field massive stars 
can quantify limitations on the minimum stellar group- 
ings needed for massive star formation. Such a study, 
which we present in this work, provides direct observa- 
tional constraints for the two competing theories of star 
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formation. 

The core accretion model suggests that stars of all 
masses form by a fragmentation process in molecular 
clouds, where cores collapsing due to self-gravity rep- 
resent the mass available to form an individual star or 
stellar multiple system (e.g. Shu et al. 1987). In this 
model, massive stars must necessarily form from mas- 
sive cores; however, it is unclear how such massive cores 
(up to hundreds of Jeans masses) can collapse without 
further fragmentation. Analytic models by Krumholz 
k, McKee (2008) suggest that sufficiently dense clouds 
having surface densities > 1 g/cm~ 2 will trap stellar 
and accretion radiation that heats clouds and prevents 
further fragmentation. Additionally, 3-D hydrodynamic 
simulations by Krumholz et al. (2009) reveal that self- 
shielding occurs along filaments resulting from gravita- 
tional and Rayleigh- Taylor instabilities, thereby channel- 
ing gas onto massive stars despite radiation pressures 
that dominate gravitational forces. These simulations 
result in the formation of a high-mass star or multiple 
system with a small companion population of low-mass 
stars. Similarly, Spaans & Silk (2000) show that the 
star formation properties of a gravitationally collapsing 
molecular cloud are highly dependent upon the equation 
of state of that cloud. For a cloud polytropic equation 
of state given by P oc p 1 where P is the thermal pres- 
sure and p is the gas density, they concluded that 7 > 1 
yields a peaked stellar initial mass function (IMF) rather 
than a power-law distribution. Li et al. (2003) conduct 
further simulations and find that molecular clouds with 
7 > 1 will most likely result in the formation of massive, 
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isolated stars. 

In contrast, the competitive accretion model suggests 
that fragmentation only produces low-mass stars, with 
high-mass stars formed by winning a competition for the 
remaining gas (e.g. Zinnecker 1982). In this scenario, the 
mass of a star is highly dependent upon the star-forming 
environment, with high-mass stars preferentially located 
at the bottom of the gravitational potential where the 
majority of a cluster gas reservoir gets funneled (Bon- 
nell et al. 2001). This model of massive star formation 
requires that massive stars form in a clustered environ- 
ment, with an explicit relation between the mass of a 
cluster (M c i) and the mass of the most massive star in 
the cluster (r7i max ) given by M c \ oc m^ x (Bonnell et al. 
2004). Thus, competitive accretion forms massive stars 
along with a fully populated cluster of lower mass com- 
panion stars (Bonnell et al. 2007). 

One of the primary differences between the observa- 
tional predictions of these models is in the formation of 
high-mass stars in low M c \ environments. Competitive 
accretion argues that formation of a high-mass star in a 
low-mass cluster is extremely difficult, while core accre- 
tion places no formal constraint on cluster mass. The 
competitive accretion model implies that the IMF is not 
a universal property of star formation, but instead tends 
to limit m mal for a given M c \. However, the IMF has 
been robustly verified for a wide range of star-forming 
environments, leading many to argue that m up , the up- 
per stellar mass limit, is a universal property of star for- 
mation, regardless of environment (see Elmegreen 2000, 
2006, 2008). In the case of a universal IMF, the relation- 
ship between m max and M c \ is determined by the statisti- 
cal mean (Oey & Clarke 2005). Not everyone agrees with 
the universality of the IMF; for example, WK06 argue for 
the existence of a deterministic m max -M c i relation using 
analytic models. They further their argument by aggre- 
gating a data sample of Galactic clusters from which they 
find a strong correlation between m max and M c \, closely 
following the derived m max -M c i relation in Weidner & 
Kroupa (2004). Weidner et al. (2010a), who update and 
greatly expand the observational data set, conclude us- 
ing a variety of statistical tests, that it is highly unlikely 
that the sample of Galactic clusters is generated from 
random sampling of a universal IMF. 

However, Selman & Melnick (2008), using the same 
data from WK06, argue that the correlation of m max 
with M c \ may be caused by the quick dispersal of clusters 
dominated by a single massive star due to gravitational 
instabilities. Since these objects would no longer be iden- 
tifiable as clusters, such a dispersal effect would bias the 
WK06 cluster sample against clusters that formed with a 
flatter-than-Salpeter mass function, leaving behind only 
those clusters that follow a more standard Salpcter mass 
function. Maschberger & Clarke (2008) complement the 
WK06 data set with a sample of very small clusters from 
Testi et al. (1997) and find that the resultant ensemble 
of clusters does not significantly deviate from the expec- 
tations of a universal stellar IMF, when examining the 
correlation between the number of stars in a cluster (iV«) 
and m max . They argue that analyses using iV* instead 
of M c \ are more reliable since AT* is a directly observable 
quantity, while M c \ must be inferred. They caution that 
observational and sample selection effects can greatly in- 



fluence the correlation of m max with M c \ or TV* and that 
much more observational data is needed to reach a con- 
clusion. 

A m max -M c i relation, if it exists, has broad implica- 
tions for cumulative stellar populations of galaxies. The 
power-law form of the cluster mass function is robust, 
similar to the IMF, with empirical derivations from a 
wide range of Galactic objects and environments, gener- 
ally consistent with a power-law slope of -1.7 to -2.3 
(§3.2). Coupling the cluster mass function, which is 
highly weighted towards low-mass clusters, with a de- 
terministic m max -M c i relation can have a large effect on 
the integrated galactic initial mass function (IGIMF) for 
stars. The primary consequences include a decrease in 
the expected number of OB stars within galaxies and 
an overall steepening of the IGIMF for the compos- 
ite stellar population of a galaxy (Kroupa & Weidner 
2003; Weidner & Kroupa 2005). A steepened IGIMF 
appears to successfully reproduce a variety of poorly- 
understood observationally-derived relations, including 
the dwarf galaxy mass-metallicity relation (Koppen et 
al. 2007), global correlations between Ha to UV flux 
ratios and galaxy mass (Hoversten & Glazebrook 2008; 
Lee et al. 2009; Meurer et al. 2009), and sharp radial 
surface brightness truncations in Ha compared to more 
extended-UV emission in the outer disks of nearby galax- 
ies (Thilker et al. 2007; Pflamm-Altenburg et al. 2009). 
Such observational relations appear to arise naturally 
from clustered star formation and the m max -M c ] relation 
implicit to a steepened IGIMF. Maschberger et al. (2010) 
found evidence of a steepened IGIMF in the competitive 
accretion simulations of Bonnell et al. (2003, 2008), link- 
ing these two theories under their common assumption of 
a m max -M c i relation. Considering the far-reaching impli- 
cations of a steepened IGIMF, it is of utmost importance 
to examine its validity using observational constraints of 
isolated O stars. 

One observational method to test the assertion of a 
w max -M c i relation is to look for isolated, massive star 
formation. Field O stars are abundant in the literature 
(e.g., Massey et al. 1995) and may account for 25-30% 
of the O star population in a galaxy (Oey, et al. 2004). 
While many of these stars are likely to be runaway stars 
from clusters, the remainder of field stars with no ev- 
idence of companions would be difficult to incorporate 
into the m max -M c i relation proposed by WK06 and in- 
herent to the theory of competitive accretion. In a study 
of Galactic field O stars, de Wit et al. (2004, 2005) find 
that 4 ± 2% of all Galactic O stars appear to have formed 
in isolation, without the presence of a nearby cluster or 
evidence of a large space velocity indicative of a runaway 
star. This value is in agreement with the 5% of isolated 
O stars (defined as O stars without any companion O or 
B stars) found from Monte Carlo simulations of clusters 
(Parker & Goodwin 2007). 

In this paper, we examine the stellar environment 
around field O stars to probe the limiting cases where 
O stars form in the sparsest stellar groups, or in near 
isolation. 

2. OBSERVATIONS 
2.1. HST Imaging Observations 
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We target field O stars in the SMC for this study be- 
cause this nearby galaxy offers a view unobscured by gas 
and dust, allowing clear identification of the field mas- 
sive stars and any low-mass companions. Our targets are 
taken from the work of Oey et al. (2004), who applied 
a friends-of-friends algorithm to photometrically iden- 
tify OB star candidates, thereby identifying clusters and 
field stars in this sample. For this study, all of the tar- 
gets were spectroscopically verified as O or early B stars, 
and all appeared isolated in ground-based imaging. In 
a pilot SNAP program, we exploit the 0.05"/px spatial 
resolution of the Advanced Camera for Surveys (ACS) 
Wide-Field Camera aboard the Hubble Space Telescope 
(HST) to search for low-mass stars associated with the 
target OB stars. Unfortunately, Cycle 14 had an unusu- 
ally low SNAP return, and we obtained observations of 
only eight targets. Table Q] lists our sample. Column 1 
gives the star's ID from the catalog of Azzopardi & Vi- 
gneau (1975); the star smc-16 was catalogued by Massey 
et al. (1995). Columns 2, 3, and 4 list the right ascen- 
sion, declination, and V magnitude, respectively, taken 
from Massey (2002). Column 5 gives the spectral types, 
in some cases derived from our observations described 
below in H2.3I Column 6 gives our mass estimate derived 
from the spectral type as described below, in ^2.41 Col- 
umn 7 gives our measured heliocentric radial velocities 
(see §2.3). 

We obtained exposures of 6 seconds in the F555W band 
and 18 seconds in the F814W band. Figure Q] shows the 
F814W images of each object, with a circle of radius 
one parsec (3.4 arcsec) for reference, adopting an SMC 
distance of 60 kpc (Harries et al. 2003). The F555W 
exposures are complete down to 21st magnitude while 
the F814W are complete down to 22nd magnitude. 

For each HST field, we use the IRAF DAOPHOT pack- 
age to identify stars and obtain their photometry, using 
a combination of aperture photometry and PSF fitting. 
The aperture photometry is more reliable due to an un- 
dersampled PSF; however, PSF fitting was unavoidable 
in the case of close companions. The aperture photome- 
try was done with an aperture radius of six pixels, while 
PSF fitting was done with a two-pixel FWHM and cor- 
rected to match the six-pixel aperture. For the F555W 
images, 3-5% of stars require PSF fitting, while in the 
F814W band, 10-15% require PSF fitting. PSF fitting 
did not reveal close companions for any of the target 
stars; however, the PSF subtraction of the target stars 
was very non-uniform, so companions cannot be entirely 
ruled out. 

We used two separate methods to identify possible 
companion stars associated with the targets: (1) an anal- 
ysis of the stellar density surrounding the OB star and 
(2) a friends-of-friends algorithm. For the first method, 
we computed the stellar surface density of the field as 
a function of radius from the target OB star. We per- 
formed this analysis using the F814W exposures, since 
they probe to a fainter magnitude than the F555W im- 
ages. To measure the average stellar density of the field, 
we used an annulus centered on the OB star, with an 
inner radius ranging from 5" to 10", and an outer ra- 
dius 10" beyond the inner radius; the annuli were posi- 
tioned to avoid any obvious stellar clustering. We define 
a density enhancement to occur when the observed stel- 
lar density is higher than the average stellar density of 



the field, including the statistical uncertainty. The prob- 
ability that the observed stellar distribution matches the 
expected Poisson distribution, for a given background 
stellar density, is: 

f(k,x) =x k *e- x /k\ (1) 

where x is the expected star count within a given radius 
and k is the observed star count within that radius. A 
more useful value, which we designate the "field proba- 
bility" P(f), is the likelihood that the population follows 
the stellar density of the field: 

k 

p(/) = i-£/(m) ( 2 ) 

o 

The right term in equation [2] yields the probability of 
observing more than k stars within a certain radius, and 
so P(f) represents the probability of obtaining the back- 
ground field. A smaller P(f) therefore indicates an in- 
creased likelihood of clustering. 

To examine the stellar environment near each target 
star, we plot the cumulative stellar density as a function 
of radius from the target OB star in Figure [2] The hor- 
izontal lines show the measured background density of 
each field. We were unable to detect stars within the 
wings of the target star, which typically extended to a 
radius of 0.2 pc. Thus, the actual cumulative stellar den- 
sities may be higher than those observed. Table [2] lists 
the angular and physical radius R c \ at which P{f) is 
minimized in columns 2 and 3, respectively. The value 
of P(f) is listed in column 4 for targets showing a density 
enhancement. 

We also searched for density enhancements with a 
friends-of-friends algorithm applied to our F814W im- 
ages. This method defines group members to be all 
the stars within a fixed clustering length I of another 
member of the same group. Following Battinelli (1991), 
we adopted a value for I that maximizes the number of 
groups having at least three stars. The distributions of 
clusters vs. I for each field are approximated well by 
normal distributions, and so wc used gaussian functions 
to estimate I. Figure [3] shows a representative example. 
The average I for these observations is 0.48 ± 0.05 par- 
sees. Table [5] lists the clustering length of each field in 
column 5. 

2.2. Minimal O Star Groups 

As seen from the data in Figure [5] and Table [5J there 
are three stars that show robust evidence of small, asso- 
ciated stellar groups: AzV 67, AzV 106, and AzV 302. 
The field probabilities, P(f), for these targets range from 

0. 01% to 11.4%. The friends-of-friends algorithm con- 
firms the existence of stars within the clustering length 
I in each of these fields. To examine the sensitivity to 

1, we also ran the friends-of-friends code with values of 
/ ± 0.05 pc, where this variation is equal to the standard 
deviation of I in the sample (see above) . For the smaller 
values, this yielded companions only for AzV 302; while 
the larger value resulted in associated stars for all tar- 
gets in our sample. Thus, we are confident that the fit- 
ted peak values for I are appropriate, and they support 
the identification of groups found by identifying density 
enhancements. We will refer to these sparse groups as 
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Table 1 

SMC Field OB Stars 



Field OB Star 


RA a 


Dec a 


ya 


Spectral Type b 


Mass (Mq) 


RV b (km s" 1 ) 


smcl6 


01:00:43.94 


-72:26:04.9 


14.38 


09 V 


23 ±2 


121 ± 21 


AzV 58 


00:49:57.84 


-72:51:54.4 


14.29 


B0.5 III 


22 ±2 


146 ± 11 


AzV 67 


00:50:11.13 


-72:32:34.8 


13.64 


08 V 


37 ±3 


159 ± 13 


AzV 106 


00:51:43.36 


-72:37:24.9 


14.18 


Bl II 


18 ±1 


150 ± 12 


AzV 186 


00:57:26.99 


-72:33:13.3 


13.98 


08 III((f)) c 


33 ±3 


159 ± 10 


AzV 223 


00:59:13.41 


-72:39:02.2 


13.66 


09 II d 


32 ±2 


189 ± 7° 


AzV 226 


00:59:20.69 


-72:17:10.3 


14.24 


07 III((f)) c 


35 ±3 


146 dz 21 


AzV 302 


01:02:19.01 


-72:22:04.4 


14.20 


08.5 V c 


27 ±2 


161 ± 11° 



a From Massey (2002). 

b Observed with the IMACS multislit spectrograph on the 6.5m Magellan /Baade telescope, unless otherwise stated. 
c From Massey (2009). 
d From Evans (2006). 

° Observed with the MIKE echelle spectrograph on the 6.5m Magellan/Clay telescope. 



Table 2 

Cluster and Stellar Population 



Field 



Angular Size (") 



Rd (pc) 



P(f) 



I (pc) 



Field IMF 



smcl6 
AzV 58 
AzV 67 
AzV 106 
AzV 186 
AzV 223 
AzV 226 
AzV 302 



10.3 
6.9 



1.8 



1.5 
1.0 



0.7 



0.114 
0.047 



0.0001 



0.52 
0.43 
0.50 
0.43 
0.46 
0.47 
0.48 
0.58 



F : 

r : 
r : 
r : 
r : 
r : 
r : 
r : 



-0.9 ±0.5 
-1.2± 1.1 
-1.3 ±0.7 
-1.6 ± 1.3 
-0.9 ± 1.2 
-1.1 ±0.5 
-1.0± 1.0 
-1.0 ±0.7 



"minimal O star groups". AzV 186 is the only field to 
appear clustered using the friends-of-friends algorithm, 
but not the density enhancement algorithm. We iden- 
tify only those fields that appear clustered using both 
algorithms, as minimal O star groups. 

2.3. Isolated Field Stars 

The remaining four targets (AzV 58, AzV 223, AzV 
226, smcl6) show no evidence of associated stars using 
either of the methods above. These are candidates for 
massive stars that formed in complete isolation. How- 
ever, a substantial fraction of field O stars likely did not 
originate in the field at all, but rather are runaway stars 
from clusters (e.g., Blaauw 1961; Gies 1987; Hoogerwerf 
et al. 2001). Thus, we examine the likelihood that our 
isolated OB stars are runaways. 

We are performing a complete spectroscopic survey of 
SMC field OB stars that is now underway, using pri- 
marily the IMACS multi-slit spectrograph on the 6.5- 
m Magellan/Baade telescope at Las Campanas Obser- 
vatory (Lamb et al, in preparation). These data yield 
the stellar radial velocities and spectral classifications, 
and we will ultimately estimate the runaway fraction for 
these massive field stars. In the course of this survey, we 
obtained spectra of the target field OB stars in our HST 
imaging sample. Table Q] gives our measured heliocen- 
tric radial velocities and spectral classifications. Several 
stars were previously classified, and we either confirmed 
or revised the spectral types, as shown. 

We identify two of these stars, smcl6 and AzV 223, as 
runaways, defined as stars having > 30 km s _1 difference 
(de Wit et al. 2005) from the SMC systemic velocity of 



155 km s" 1 (Staveley-Smith et al. 1997). We note that 
Stanimirovic et al. (2004) found the SMC to have a ve- 
locity gradient; however, the positions of our targets do 
lie between the 150 km s _1 and 160 km s _1 contours as 
plotted in Figure 3 of Stanimirovic et al. (2004) . The two 
runaway stars will be removed from analysis in all subse- 
quent sections, since they are no longer located in the re- 
gion of their formation. This leaves us with two, isolated, 
apparently non-runaway OB stars, AzV 58 and AzV 226. 
Assuming an isotropic distribution of runaway velocities, 
we still expect to miss a number of transverse runaways. 
The missed fraction depends on the velocity distribution 
of runaway stars; however, we estimate that typical ejec- 
tion velocities of 60 km s _1 or 120 km s^ 1 would cause 
us to miss two or one transverse runaway (s), respectively. 
Therefore, transverse runaways may account for both our 
remaining isolated stars. 

It is also informative to investigate the interstellar gas 
around these field OB stars (Figure 2} . We examined the 
ionized gas around our target stars, using Ha data from 
the Magellanic Cloud Emission Line Survey (MCELS; 
Smith et al. 2000). Since gas is a necessary component of 
star formation, the presence of gas can help to constrain 
which stars may still be in the region of their formation. 
As a control group, all three of our minimal O star groups 
show Ha emission within 2.0' (35 pc), consistent with a 
physical association within these groups as sparse, young 
clusters. We also find that both confirmed runaway stars 
are far removed from any Ha emission. For the remaining 
isolated targets, the MCELS data show that AzV 58 and 
AzV 226 are located within Hn regions in the line of sight 
(Figure S]), while AzV 186 is far from any Hn regions. 
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Figure 1. HST F814W images of each target star, with the star's ID number from the Azzopardi & Vigneau (1975) catalog listed in the 
upper left of each panel. The circle corresponds to a radius of 1 parscc. The top row contains fields showing a stellar density enhancement. 
The middle and bottom rows contains fields with no density enhancement. The bottom row contains stars found to be runaways. 



These results suggest that AzV 58 and AzV 226 may still 
be in the region of their formation and thus they remain 
candidates for isolated massive star formation. 

2.4. The IMF 

We construct color-magnitude diagrams (CMD's) from 
our ACS photometry, and use these to evaluate the IMF 
for the minimal O star groups as well as the field stars. 
We converted the photometric results from the F555W 
and F814W bands to Johnson Fand Cousins /bands, re- 
spectively, following the synthetic transfer equations pro- 
vided in Sirianni et al. (2005). These V and I band mag- 



nitudes were then extinction corrected using the SMC 
extinction maps provided by the Magellanic Cloud Pho- 
tometric Survey (MCPS; Zaritsky, et al. 2002). MCPS 
provides two sets of stars for estimating extinction, the 
"hot" stars (12, 000 K < T eff < 45, 000 K) and "cool" 
stars (5, 500 K < T eff < 6, 500 K). We chose to average 
the extinction values of "cool" stars within V of the tar- 
get star for the extinction calculation, as these objects 
better trace the sparse field population of our observa- 
tions. The "hot" stars tend to trace active, clustered star 
formation with high gas content, yielding extinction val- 
ues up to half a magnitude higher than those typically 
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Field 67, 08 V 



4 r 



3 r 





0.0 



1.0 1.5 
Radius (pc) 

Field 186, 08 li(f) 




0.5 



1.0 1.5 
Radius (pc) 



2.0 




1.0 1.5 
Radius (pc) 




1.0 1.5 
Radius (pc) 



Figure 2. Cumulative stellar density as a function of radius from each target star. Each field is titled with the target star's ID and 
spectral type. The panels are shown in the same sequence as in Figure[l] with the first three stars showing a stellar density enhancement. 
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Field 67 




Clustering Length (pc) 

Figure 3. Number of clusters vs. clustering length, I, with a fitted 
Gaussian (dashed line) for a representative field, AzV 67. The bar 
shows the extent between one standard deviation below and above 
the peak value, and indicates the three values used for I for each 
field. 



found in the "cool" stars. 

We then calculated absolute magnitudes using a dis- 
tance modulus of 18.9 (Harries et al. 2003). Compar- 
isons of our photometry with MCPS photometry typ- 
ically show differences less than 0.1 mag, with /-band 
matching more closely than V. On average, our pho- 
tometry appears fainter than that of the MCPS in both 
V and I bands. For stars fainter than 16th magnitude, 
Zaritsky et al. (2002) find that stellar crowding may play 
a factor in spuriously brightening the MCPS photometry 
when compared with the OGLE catalog (Udalski et al. 
1998), which they deem to be superior in this regime. 
Since nearly all stars in our fields are fainter than 16th 
magnitude, the difference between our photometry and 
MCPS photometry is likely due to these crowding effects. 

We used Geneva stellar evolutionary tracks to extract 
mass estimates from our photometry. These tracks are 
based on the ATLAS9 no-overshoot models at SMC 
metallicity, calculated by Charbonnel et al. (1993) and 
have been converted by Girardi et al. (2002) to Johnson- 
Cousins UB VRI photometry. The stellar masses are in- 
ferred by identifying the two evolutionary tracks between 
which a star falls, on the V vs V — I CMD. For the ex- 
tremely blue, O and B target stars, the V — I colors are 
degenerate, so effective temperature, T e g , is used instead 
ofV — I color. T e ff of our OB target stars are based on 
our spectral classifications given in Table [TJ using the 
conversions of spectral type to effective temperature for 
the SMC by Massey et al. (2005) for the O stars, and 
Crowther (1997) for the B stars. 

Following the formalism of Scalo (1986) and analysis 
done by Massey (1995), we write the slope of the IMF 
as: 

r = d log £(log m) jd log to (3) 

where £(log to) is the mass function in units of stars born 
per logarithmic stellar mass to (Mq) per unit area (kpc 2 ) 
per unit time (Myr). This corresponds to a power- law 
mass spectrum given by: 

n(m) dm oc m? 'dm (4) 



where n(m) dm is the number of stars per unit mass bin 
and the power law index 7 = T — 1. In this formalism, a 
Salpeter mass function has a slope T = —1.35. In order 
to derive £(log to), we counted the stars in each mass bin, 
corrected for the size of the mass bin by normalizing to 
one dex in mass, and divided by the area covered by the 
observations. For computing the IMF of the background 
field, we divided by the average age of stars within each 
mass bin (Charbonnel et al. 1993) to account for dif- 
ferences in stellar lifetimes as a function of mass. Ages 
were calculated as an average between the lower and up- 
per mass in each bin, weighted by the IMF. This allowed 
us to measure an IMF under the assumption of continu- 
ous star formation, rather than obtaining the present-day 
mass function of the field. 

The CMD and IMF for a representative, full field, that 
of AzV 67, are shown in Figure [5] On the CMD plot, we 
draw SMC-metallicity evolutionary tracks (Charbonnel 
et al. 1993; Girardi et al. 2002) to show the stellar mass 
ranges of the field population. On the IMF plot, the er- 
ror bars represent the Poisson uncertainty for each mass 
bin. The detection limit is V = 22 in the F555W image, 
which corresponds to a mass of 1.5M Q . The observations 
are incomplete up to 2M Q . On the IMF plot, points cor- 
responding to stellar mass bins below 2M Q are excluded 
from the linear fitted line, which has been made weight- 
ing the data inversely by the errors. The IMF slopes for 
the field population in each full, 202" x 202" ACS frame 
are given in Table [21 column 6. Accounting for the stel- 
lar age correction, each field exhibits an IMF consistent 
with a Salpeter IMF, within the uncertainty. 

Due to the small number of cluster members identified 
in both the F814W and F555W exposures, it is not use- 
ful to derive IMFs for the three minimal O star groups 
individually. Instead, we created a composite CMD of 
all members from the three groups. Because of field 
star contamination, only a subset of the stars compris- 
ing these density enhancements are physically associated; 
the background field star density implies a contamina- 
tion of 1 - 2 stars per target frame. In deriving the IMF 
of these minimal O star groups, we exclude red giants, 
since they are certainly field stars unassociated with re- 
cent star formation. Some main sequence stars may also 
be field star contaminants, but they are indistinguishable 
from true cluster members. Table [3] provides photometry 
of the main sequence companion stars present in both V 
and / images, above the completeness limit of V = 21 
and I = 22 magnitudes. Columns 1 and 2 list the right 
ascension and declination, while Columns 3 and 4 list the 
V and I magnitudes. 

Figure [6)3 shows the CMD for all companions found 
using the density enhancement analysis. The CMD for 
companions found using the friends-of-friends algorithm 
is identical to that in Figure [Bh above our completeness 
limit of 2Mq. We plot the composite IMF of these min- 
imal O star groups in Figure O), excluding all mass bins 
below our completeness limit. We find one companion 
star in each of the mass bins, 2M < to < 2.5M Q , 
2.5M Q < to < 3M Q , and 3M < m < 4M Q . Upon 
generating IMFs of the companions found in the friends- 
of-friends analysis with clustering lengths I i la, we 
find that changing the exact prescription for determining 
companionship has little effect on the measured IMF for 
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Figure 4. Ha images from the Magellanic Cloud Emission Line Survey (MCELS), tracing the ionized gas in the region around each target 
star. The PSF of these observations is similar to the 1 parsec circles in Figure [l] and the field of view is 9' X 9'. The panels are sequenced 
as in Figure[T] with the top row showing fields with a stellar density enhancement. The target stars are identified. 



this population of minimal star groups. 

The IMF slope for the composite population in the 
minimal star groups varies between T = 0.1 ± 1.0 to 
— 0.2±0.9 in the preceding analysis, which is not far from 
a Salpeter slope of T = —1.35 within the uncertainties. 
However, the true slope may be even larger, since we 
do not correct for contamination by main sequence field 
stars. Due to the presence of very massive stars and a 
small number of low-mass companions, the IMF is es- 
sentially predetermined to be unusually shallow, but we 
note that it does not necessarily represent a significant 
variation from the Salpeter value, statistically, given the 



small numbers of stars. 

3. MONTE CARLO SIMULATIONS 

This regime of sparse star formation is strongly dom- 
inated by stochastic effects and, as described in §1, it 
offers an important discriminant between star formation 
theories. A fundamental question is whether the stellar 
IMF in clusters is largely independent of parent cloud 
mass, and determined by only the most local physics (e.g. 
Maschberger & Clarke 2008; Krumholz et al. 2010). If 
so, the IMF in the lowest-mass clusters remains identical 
to that seen in higher-mass clusters, and should be de- 
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67 Field Star CMD 67 Field IMF 




V-l log(m [M ]) 

Figure 5. (a) Color-magnitude diagram and (6) field IMF for a representative field, AzV 67. Each of our observed fields is consistent 
with the Salpctcr IMF within the uncertainty. The error bars represent the Poisson uncertainties. Mass bins below our completeness limit 
of 2Mq are excluded from the fitted line. 



Table 3 

Photometry of Companion Stars 



O Star Group 


RA 


Dec 


V 


I 


AzV 67 


00:50:10.39 


-72:32:30.9 


20.78 ±0.15 


20.77 ±0.06 


AzV 67 


00:50:11.77 


-72:32:30.0 


19.85 ±0.07 


19.79 ± 0.03 


AzV 67 


00:50:11.95 


-72:32:32.3 


20.51 ±0.13 


20.67 ±0.06 


AzV 106 


00:51:43.45 


-72:37:27.3 


19.51 ±0.05 


19.58 ±0.02 


AzV 106 


00:51:44.02 


-72:37:23.9 


21.20 ±0.20 


21.15 ±0.08 


AzV 302 


01:02:18.74 


-72:22:04.7 


19.76 ± 0.08 


19.52 ±0.03 


AzV 302 


01:02:18.81 


-72:22:01.1 


20.84 ±0.16 


20.68 ±0.06 


AzV 302 


01:02:18.83 


-72:22:01.7 


17.97 ±0.02 


17.98 ±0.01 


AzV 302 


01:02:19.40 


-72:22:02.2 


20.81 ±0.16 


20.46 ±0.04 



Companion CMD Composite IMF 




0.0 0.5 1.0 1.5 0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 

V-l log(m [M ]) 



Figure 6. Shown in panel (a) is the CMD of the minimal O star groups with membership determined by the stellar density analysis. 
Panel (b) shows the IMF of the minimal O star groups including companions down to our completeness limit of 2Mq for the F555W band. 
We plot a solid line, fit to the data, with T = 0.1 ± 1.0. For reference, a dashed line with V = —1.35, representing a Salpeter slope, is shown. 
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scribed by the simple random drawing of discrete stars 
from an IMF represented as an ordinary probability den- 
sity function. On the other hand, a scenario that is at 
least as plausible is that the IMF is driven by, and limited 
by, the mass of the parent molecular cloud (e.g. Bonnell 
et al. 2004; Weidner & Kroupa 2004). If so, then the 
stochastic effects seen in the lowest-mass clusters will be 
linked to the parent cloud masses. The objects in our 
study offer a unique opportunity to explore this extreme 
parameter space. 

3.1. Analytic Probability 

We can quantify the likelihood that these minimal O 
star groups conform to a standard but under-sampled 
stellar IMF. For a cluster of stars, the probability 
that all stars are below mass m max is given by: 



-P(m max , AT») 



(m) dm 



(5) 



where 0(m) is the IMF and m m ; n is its lower mass limit, 
which we assume to be constant. We adopt a Kroupa 
IMF (Kroupa 2001), whose form is: 



4>{m)dm oc 



m 



r 1.3 
-2.35 



dm , O.O8M < m < 0.5M o 
dm , O.5M <m< 15OM 



(6) 



For the minimal O star groups, we estimated as fol- 
lows. We count the candidate member stars, identified 
as described in §§2.1 - 2.2, and correct for the expected 
number of field stars contaminating the cluster line of 
sight. The contamination is determined by the stellar 
density of the background held as calculated in i )2.1l and 
the angular size of the cluster (Table [2]). We further cor- 
rect the observed numbers to account for stars below our 
completeness limit of 1.5M Q in the F814W images to get 
a hnal estimate over the full mass range of the IMF. 
Table2]hsts star counts for each held as follows. Column 
1 shows the star ID; columns 2 and 3 show the total num- 
ber of stars observed within the cluster raduis and the 
subset of those estimated to be field stars, respectively; 
column 4 lists the resulting number of cluster members 
above the detection threshold; and column 5 lists the 
inferred N*, integrating over the full stellar mass range 
from m\ Q = 0.08M Q . We follow a similar process to esti- 
mate cluster mass by summing the masses of all stars in 
the cluster and again integrating below our completeness 
limit down to m\ — 0.08Mq. The cluster mass estimates 
are listed in column 6 and will be utilized for analysis in 
§5- 

We then use these A* values and the m max values from 
Tabled] to calculate P(m max , N*) (equation [5]) . These 
values are given in column 7 of Table |H For the minimal 
O star groups, we find that the likelihoods of these clus- 
ters containing stars as massive as those observed range 
from 11 to 20%. In the case of those isolated O stars with 
no observed companions, and assuming that all cluster 
members are below our detection threshold, the likeli- 
hoods that the observed stars formed with no compan- 
ions above 1.5Af© ranges from upper limits of <1% to 
<3% (Table H]). These likelihoods, while low, are not 
exceedingly so, suggesting that the occurrence of these 
minimal O star groups, even those with no stars above 
the detection threshold, is not especially unlikely given 



the assumed parameters. Our target selection was based 
on apparent isolation in ground-based imaging, and so we 
expect our objects to fall in this low-probability regime. 
As a reminder, we note that these probabilities are based 
on maximized numbers of stars below the detection limit, 
in all cases. 

3.2. Numerical Simulations 

We now devise Monte Carlo simulations that generate 
clusters and cluster members to explore the frequency of 
the observed minimal O star groups in the context of a 
stellar clustering law or cluster mass function. 

We simulate the scenario for which the IMF is com- 
pletely independent of cluster mass. We generate clus- 
ters using a stellar clustering law having a default (3 = — 2 
power law slope (hereafter the "iV* simulations" ) : 



n(N m ) dN* oc N? dN* 



(7) 



where n(JV») dN* is the number of clusters in the range 
to + dN*. We set single stars, N* — 1, as our 
minimum "cluster" and set an upper limit for clusters 
at N„ — 10 6 . Our use of the -2 power law index is mo- 
tivated by observations of a wide variety of stellar pop- 
ulations including young, massive clusters (e.g., Hunter 
et al. 2003; Zhang & Fall 1999; Fall et al. 2009), super 
star clusters (Meurer et al. 1995), globular clusters (e.g., 
Harris & Pudritz 1994), and HII regions (Oey & Clarke 
1998). Oey et al. (2004) found that this power law ap- 
plies smoothly down to individual field OB stars, and we 
assume that the same clustering law holds true for all 
stellar masses, down to our lower mass limit of 0.08-Mq. 

Once clusters are generated following the A* clustering 
law (equation!?]) , each cluster is randomly populated with 
stars using the IMF given by equation [6] The IMF, in- 
cluding its stellar mass limits, is constant for all clusters, 
thereby allowing true, isolated O stars to be generated 
in A* — 1 "clusters" . For each AT* simulation, up to 10 7 
clusters are generated to ensure that we create enough 
clusters to fully populate the cluster parameter space. 

We also carry out simulations that generate clusters by 
total cluster mass M c \, according to essentially the same 
power-law distribution (hereafter the U M C \ simulations" ) : 



n(M cl ) dM cl oc dM ci 



(8) 



where n{M c \) is the number of clusters in the range M c \ 
to M c \ + dM c \. We adopt an upper limit to the cluster 
mass function of 10 5 M Q . 

Most simulations of clusters heavily favor this method 
of populating clusters by mass (e.g. WK06; Parker & 
Goodwin 2007; Haas & Anders 2010). However, the ex- 
act prescription for populating stars up to the target M c \ 
varies. For our simulations, we follow Parker & Goodwin 
(2007) by populating the cluster with stars, randomly 
sampled from the IMF, until the cluster contains at least 
98% of its target mass in stars. At this point, if the last 
star generated pushes the cluster mass past 105% of the 
target cluster mass, the entire cluster is discarded. The 
process of populating the cluster with stars is repeated 
until the total stellar mass falls within 98% to 105% of 
the target cluster mass. If these conditions are not im- 
posed, the high-mass stars are too often generated as the 
last star in the cluster, and the cluster mass function is 
not well preserved in the final sample of clusters. An 
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Table 4 

Observed and Estimated Cluster Membership 



Star 


Observed 


Field 


Detected 


Estimated 


Estimated 


P(m max , N t ) 




Stars 


Stars 


Members 


TV* 


M cl (M ) 


AzV 58 


1 





1 


<19 


<43 


<0.03 


AzV 67 


24 


15 


9 


171 


103 


0.11 


AzV 106 


17 


9 


8 


152 


72 


0.20 


AzV 186 


1 





1 


<19 


<72 


<0.01 


AzV 226 


1 





1 


<19 


<77 


<0.01 


AzV 302 


13 


3 


10 


190 


81 


0.19 



investigation of the effects on cluster population using 
different algorithms can be found in WK06 and Haas & 
Anders (2010). 

We also perform variations of our simulations using 
/3 = —1.8 and -2.3 power law slopes. These are val- 
ues typically observed as empirical variation for real 
systems. For example, observations of giant molecular 
clouds (GMCs) and self-gravitating clumps within GMCs 
show power-law mass distributions with slopes /? ~ —1.7 
(see Bertoldi & McKee 1992; Rosolowsky 2005). Obser- 
vations of star clusters in the Large Magellanic Cloud 
show measured values of /? ~ —2.0 to -2.3 (Hunter et al. 

1993) and f3 1.8 (Chandar et al. 2010). Allowing /3 

to vary in our simulations allows us to explore its effect 
on our results. 

3.3. The Cluster Lower-Mass Limit 

Our cluster simulations are designed to probe the limit- 
ing case of sparse O star formation. A critical parameter 
in this regime is the lower limit placed upon cluster mass 
or membership number. 

With respect to quantized star counts, observational 
evidence for a -2 power law relation extending to A* = 1 
clusters are limited to samples of O and OB stars, prob- 
ing stellar masses > 2OM and > 10M Q (Oey et al. 
2004). The form of the stellar clustering law for a com- 
plete sample of stars below these masses is unknown. 
Our simulations extrapolate this -2 power law relation to 
stellar masses two orders of magnitude lower than those 
observed. However, it is possible that the clustering law 
has a turnover or cutoff in the low number regime, or the 
-2 power law relation may break down at a stellar mass 
higher than our lower mass limit of 0.08M Q . 

This low A* regime exists in an intermediate stage be- 
tween clustered formation and isolated formation, which 
Adams & Myers (2001) term "group formation". They 
estimate that the majority of star formation occurs in 
this group formation regime, with A* = 10-100. They 
also estimate a lower bound of A*= 36 for a group to 
evolve as a cluster, defined as having a relaxation time 
longer than the crossing time of the group. However, N- 
body simulations show that clusters in this "group for- 
mation" regime are likely to lose a significant portion of 
their stellar members on timescales of a few Myr (Bon- 
nell & Clarke 1999). Interactions that cause the stellar 
losses preferentially affect low-mass stars, thus making 
the initial A» a difficult parameter to estimate from ob- 
servations for clusters in this regime. This effect may 
be present in the embedded cluster catalog from Lada 
& Lada (2003), which restricts membership to clusters 
with A* > 35 and shows a turnover in the cluster mass 



function below 50M©. Since the typical stellar mass is 
~ O.5M , the turnover that they observe in the cluster 
mass function at 50M Q corresponds to A* ~ 100, right at 
the regime where the N-body simulations predicted stel- 
lar membership loss. Therefore, while the exact nature 
of the lower limit of the clustering law is still unknown, a 
truncation at a lower value of A^io in the range from ^35 
to 100 stars or M c \ from ~17.5 to 5OM is reasonable. 
This lower limit is similar to that employed in previ- 
ous cluster simulations, which typically range from 5Mq 
(Weidner & Kroupa, 2004) to 5OM (Parker & Goodwin, 
2007). 

With respect to the cluster mass function, extragalac- 
tic studies of complete samples of clusters have probed 
to M cl ~ 1O 3 M (Chandar et al. 2010), limiting our 
knowledge of the form of the mass function for complete 
samples below this mass. Lada & Lada (2003) probe 
clusters below 10 3 M Q by compiling a catalog of embed- 
ded clusters in the solar neighborhood. Their results are 
consistent with a M~^ 2 cluster mass function for clusters 
from 50M Q to 10 3 M Q . Below 50M©, they find a statis- 
tically significant turnover in the cluster mass function. 
This may indicate that the lower limit of our M c \ simu- 
lations should be truncated or the power law should turn 
over around a value 50M Q . 

In the next sections, we compare different simulations 
with observed statistics for SMC massive star clustering 
and isolated Galactic O stars. In this comparison, we in- 
clude three separate lower A* limits for the clustering law 
at A*j = 1, 40, and 100 stars, and two separate lower 
M c i limits of the cluster mass function at M c \ t \ — 20M Q 
and 50Mq. Since the typical stellar mass is ~ 0.5M Q , 
A^io = 40 is the appropriate equivalent lower limit for 
M c y = 20M Q and that for M cUo = 5OM is A*^ 
1 00. While these values are empirically motivated, they 
do also allow the formation of stars > 20M Q , as required 
for our purposes. In what follows, it is important to bear 
in mind that the cluster parameterizations are extrapo- 
lated beyond observed mass ranges. 

3.3.1. Comparison with SMC Clustering Statistics 

We compare the distribution of massive stars in the 
simulated clusters with the actual, observed distribu- 
tion of stars among SMC clusters. Oey et al. (2004) 
broke down the clustering of observed, photometrically- 
identified massive star candidates with a friends-of- 
friends algorithm. They identified two samples, the "O- 
star sample" and the "OB-star sample," corresponding 
to masses m > 20M Q and m > 10M Q , respectively. In 
our simulations, we consider only clusters containing one 
or more O or OB stars, defined by the same respective 
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mass ranges. Table [5] summarizes our findings on massive 
star clustering. The rows are divided into three sections 
which denote results from the A* simulations, the M c \ 
simulations, and the Oey et al. (2004) observations. For 
the simulations, N* t i or M c \ t \ is listed in column 1 and 
the power law slope is listed in column 2. Columns 3, 4, 
and 5 list the frequencies of having one, two, or more than 
two O stars, respectively, in the clusters; while columns 
6, 7, and 8 list these frequencies for OB stars in the clus- 
ters. 

Table [5] reveals three trends: (1) steepening the power 
law slope of the simulation results in an increasing frac- 
tion of O or OB star clusters containing a single O or 
OB star, (2) the fraction of OB star clusters with a sin- 
gle OB star is lower than the fraction of O star clusters 
with a single O star and (3) the ./V* simulations are more 
sensitive to the lower limit truncation than the M c \ sim- 
ulations. The hnal trend is due to our cluster population 
method. In the case where M c \ < m up , stars with mass 
greater than M c \ will not be allowed to form in such a 
cluster in the M c \ simulations. However, in the iV* sim- 
ulations, each cluster can form stars up to m up . 

In the observations by Oey et al. (2004), 61-65% of 
clusters having at least one massive star contain only 
a single massive star. We hnd best agreement with this 
fraction in our simulations following a -2 power law slope 
and a truncation of = 40 or 20M©. Table[5]also shows 
good agreement with (3 = —1.8 for the simulation 
having a truncation of AT, — 20. For f3 = —2.3, the 
steepness of the slope causes the OB cluster sample and 
O star cluster sample to behave differently enough from 
each other that they cannot both be in agreement with 
the observations, regardless of the lower limit. Similarly, 
at a truncation of AT, = 100 or M c \ = 50M©, the OB 
cluster sample and O star cluster sample never come into 
agreement simultaneously at the same power law slope. 

3.3.2. Comparison with Galactic Isolated O Star Fraction 

We also compare our simulations to the de Wit et 
al. (2005) result that 4±2% of all Galactic O stars are 
found in isolation. The simulations by Parker & Good- 
win (2007) are in agreement with this result. Since we 
modeled our M c \ simulations after their work, our re- 
sults should match theirs quite closely, although the ex- 
act parameters of the simulations differ slightly. Parker 
& Goodwin (2007) set the lower mass limit of the stel- 
lar IMF to 0.1 Mq , while our simulations use a value 
of 0.08M©. Also, Parker & Goodwin use a lower limit 
of M c \ = 50M Q , while we vary the limit as discussed 
above. For comparison with this study, we adopt the de 
Wit et al. (2005) definition of O stars as having mass 
> 17.5A/q. This comparison provides a good check for 
our results and also quantifies the effect of the lower limit, 
Mcijo or AT, jo, on the isolated O star fraction. 

Parker & Goodwin's definition of an isolated O star 
is twofold: first, the cluster contains no B stars with 
mass > 10M© and second, the total cluster mass M c \ < 
IOOMq. These constraints were set according to the de- 
tection limits of the de Wit et al. (2005) result and are 
intended to mimic the sensitivity of those observations. 
For the following analysis, we follow Parker & Goodwin's 
definition of O star isolation. The results of this analysis 
are shown in Table [6l which has rows divided between 
the A/ - * simulations, the M c \ simulations, and the de Wit 



et al. (2005) observations. N* t \ or M c \ t \ for the sim- 
ulations is listed in column 1, the power law slope for 
the simulations is listed in column 2, and the fraction of 
isolated O stars as a fraction of all O stars is listed in 
column 3. This table reveals trends similar to those seen 
in Table that steepening the power law slope increases 
the fraction of isolated O stars, and that the AT* simu- 
lations are more sensitive to the lower limit truncation 
than the M c \ simulations. Table [6] also shows that the 
fraction of isolated O stars matches very well between 
the N* and M c \ simulations having lower limits of A*^ 
= 40 and M c y = 20M©, respectively. 

For their simulation adopting a -2 power law slope, 
Parker & Goodwin find that 4.6% of O stars are isolated, 
in good agreement with the 4±2% result from de Wit et 
al. (2005). We find that many of our simulations match 
the de Wit et al. (2005) result, spanning the full range 
of lower limits and power law slopes tested, with the 
exception of A^jo = 1. Table [5] shows that the -2 power 
law simulations that best agree with the isolated O star 
fraction have the higher values of M c \ : \ — 50M© or N^ : \ 
= 100, whereas the simulations that best agree with the 
SMC O and OB star clustering have the lower values 
of M cUo = 2OM or A^ lo = 40 (Table EJ. Simulations 
following a -2.3 power law slope are least reconcilable 
with both sets of observations. 

3.3.3. Default Clustering Models 

In §£ !3.3.1l and l3.3.2l we find that a number of our simu- 
lations in both AT, and M c \ agree with observed statistics 
of SMC massive star clustering and Galactic isolated O 
stars. The AT, simulations are more sensitive to both the 
power law slope and N* t \ a , and many models agree with 
these two sets of observations (Tables [5] and |6|). But for 
the M c \ simulations, we find that only the f3 = —2.0 mod- 
els match with both sets of observations, albeit for differ- 
ent Mcijo lower limits. For the -1.8 and -2.3 power law 
slopes, none of the simulations can appropriately match 
both the SMC O and OB clustering simultaneously for a 
given N* t \ or Af c i,io, making a strong case for = —2.0 
models as the best choice for the power law slope. 

We therefore take the -2.0 power law as the base model 
for comparison between the M c \ and AT, simulations. 
From Table [5] we find that models with N*^\ = 40 and 
20Mq best match the SMC massive star clustering ob- 
servations, while Table [6] shows that models with lower 
limits at N* t \ — 100 and M c y = 50M© best match the 
isolated Galactic O star observations. Looking further, 
in Table [5] we see that the models with truncations at 
N* = 100 and M c \ = 50M© also match well with SMC O 
star clustering, but do not agree with SMC OB star clus- 
tering, with the fraction of single OB star clusters being 
off by ~ 10%. This difference is large compared to the 
fraction of isolated O stars for the models with N+ i — 
40 and M C \_\ Q = 20Af©, which agree with the observations 
within ~ 1 — 2%. Thus, we conclude that the simulations 
that best match both the SMC massive star clustering 
statistics and the fraction of isolated Galactic O stars are 
the A^* simulation with A»j = 40 and the M c \ simula- 
tion with Mcijo = 20A/q. We cannot rule out either AT* 
or M c \ in favor of the other and do not find evidence that 
either one is a more fundamental parameter. 

In the following sections, we compare these two simu- 
lations with our observations of minimal O star groups. 
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Table 5 

Fraction of O (OB) star clusters having a given number of O(OB) stars 



M clj io or 7V» ao 


P 


1 O star 


2 O stars 


> 2 O stars 


1 OB star 


2 OB stars 


> 2 OB stars 


N* Simulations 51 


1 


be 


0.75 


0.10 


0.151 


0.73 


0.11 


0.16 


10 


-1.8 


0.61 


0.15 


0.24 


0.53 


0.18 


0.29 


100 


-1.8 


0.53 


0.18 


0.29 


0.40 


0.20 


0.40 


1 


-2.0 


0.87 


0.07 


0.07 


0.85 


0.08 


0.08 


40 


-2.0 


0.70 


0.15 


0.15 


0.61 


0.18 


0.21 


100 


-2.0 


0.61 


0.18 


0.21 


0.46 


0.22 


0.32 


1 


-2.3 


0.96 


0.03 


0.02 


0.95 


0.04 


0.02 


40 


-2.3 


0.79 


0.12 


0.08 


0.70 


0.17 


0.13 


100 


-2.3 


0.69 


0.18 


0.14 


0.53 


0.23 


0.24 


M c \ Simulations' 1 


20M e 


-1.8 


0.56 


0.15 


0.29 


0.53 


0.16 


0.31 


5OM 


-1.8 


0.54 


0.16 


0.30 


0.44 


0.19 


0.37 


2OM 


-2.0 


0.67 


0.14 


0.19 


0.63 


0.16 


0.21 


5OM 


-2.0 


0.64 


0.16 


0.21 


0.52 


0.20 


0.28 


2OM 


-2.3 


0.78 


0.12 


0.10 


0.74 


0.14 


0.12 


5OM 


-2.3 


0.75 


0.14 


0.12 


0.61 


0.20 


0.19 


SMC Observed 15 






0.61±0.08 


0.19±0.04 


0.19±0.04 


0.65±0.04 


0.15±0.02 


0.19±0.02 



a Errors for simulated values are < 0.01 
b From Oey et al. (2004). 



Table 6 

Fraction of isolated O stars a 



M cl , lo or AT* jlo 


P 


Iso. Fraction b 


N* Simulations 


1 


-1.8 


0.099 


10 


-1.8 


0.029 


100 


-1.8 


0.008 


1 


-2.0 


0.330 


10 


-2.0 


0.083 


100 


-2.0 


0.022 


1 


-2.3 


0.728 


10 


-2.3 


0.212 


100 


-2.3 


0.057 


M c i Simulations 


2OM 


-1.8 


0.027 


5OM 


-1.8 


0.019 


2OM 


-2.0 


0.072 


5OM 


-2.0 


0.048 


2OM 


-2.3 


0.190 


5OM 


-2.3 


0.123 



Galactic Observed' 



0.04±0.02 



a Here, the definition of isolated O stars is from Parker & 
Goodwin (2007). 

b Errors for simulated values are < 0.01. 
c From de Wit et al. (2005). 

For this comparison we use the M cl 2 and N~ 2 simula- 
tions with Mcijo = 20M Q and N* t \ — 40, respectively, 
as our default models. 

4. STELLAR MASS RATIO m max , 2 /m ma x 

A simple and important parameter we can compare be- 
tween our observations and simulations is m maXj 2/TO max , 
the mass ratio of the second-most massive and most mas- 
sive stars in the cluster. This ratio is a directly observ- 
able quantity that we can measure for the minimal O 
star groups, and as such, provides a powerful parame- 
ter to use as a comparison between our observations and 
simulations. 

There are two populations in the simulation that are 



of interest: (1) clusters that contain one or more O stars 
and (2) clusters that contain just one O star. We ex- 
plore the full simulation parameter space of the mass ra- 
tio m maX! 2/m max as a function of m max in Figures [7] and 
[8l For Figures [7] and [U panel (a) includes simulated clus- 
ters with at least one O star in the cluster, while panel 
(b) includes only simulated clusters with exactly one O 
star. For these plots and all subsequent simulations, we 
adjust our definition of O stars to be m > 18M Q instead 
of our earlier 20M Q definition, allowing our lowest-mass 
target stars to be included in the parameter space when 
calculating percentile frequencies below. Since the un- 
certainty in our stellar mass estimates is 1M Q to 3M Q , 
our adjusted definition is on the order of the uncertainty. 
The simulated clusters plotted in these figures are the O 
star clusters taken from a random sample of 573 sim- 
ulated OB clusters, the same number of OB clusters 
as observed in the SMC, having at least one OB star 
(m > 10M Q ; Oey et al. 2004). In these and all subse- 
quent figures, we exclude clusters without O stars from 
the plots. The color coding in Figures [7J and [5] indicates 
the number of stars per cluster, with black indicating 
< log N* < 1, red indicating 1 < log < 2, orange in- 
dicating 2 < logiV* < 3, green indicating 3 < log < 4, 
and blue indicating 4 < log TV, < 5. Observations from 
this paper are plotted as black squares. For our appar- 
ently isolated O stars, we note that undetected compan- 
ions may exist, having individual masses up to 1.5M Q , 
and so we plot their m maX! 2/TO max as upper limits. 

Our observations lie in an interesting region of the pa- 
rameter space in both the M c \ and -/V* simulations. There 
is a strong drop-off in the population of simulated clus- 
ters having m maX! 2/m max < 0.02. This drop-off coincides 
with some of our observed isolated targets. Thus, the up- 
per limits in ^max,2/^max 

for the isolated stars approach 
the extreme lower limit of the parameter space. All of 
our observations lie within the parameter space covered 
by simulated clusters. 

To examine this quantitatively, we identify the sub- 
set of simulated clusters that correspond to the selection 



14 Lamb, J. B., et al. 







N,~ 2 , N,^40, N ^1 








NT 2 , N,>40, N =1 




1 .0 


, , , | , , , | , ., , | , , , | , , A | , , , t i i- ■ j j 

( Q ) ->>. : . ■ . 




1 .0 


- (b 


) : 




0.8 






0.8 






I 






i 








E 


0.6 




E 


0.6 




* m ■ 


1 
E 


0.4 




1 
E 


0.4 








0.2 


• . • 




0.2 








0.0 


TV- '• " '•* ". ■ ■ 




0.0 







20 40 60 80 100 120 140 20 40 60 80 100 120 140 

rn mo>! [M ] m mQ>! [M ] 

Figure 7. m, m ax,2/«imax vs m max from the N* 2 simulation with N t i D = 40 for (a) all clusters having at least one O star, and (b) all 
clusters having only a single O star. Data are color coded in logarithmic bins of N*, with black = < logN* < 1, red = 1 < logiV, < 2, 
orange = 2 < logN* < 3, green = 3 < log TV* < 4, and blue = 4 < log TV* < 5. Our observations are plotted as black squares. 



-\ N ^1 M c r 2 , N =1 



1 .0 


:(□)• 

* • * # 


i • _ • 




1 .0 






0.8 






E 


0.8 


i 




0.6 


" *" ■ 

3 V. : " " ' 




E 


0.6 


• » 

:;'«•. •. 




0.4 






E 


0.4 










E 




















0.2 








0.2 












y£ •:* • 




0.0 


1 /TV;/ 


• » 




0.0 


: 'TV; ; 




20 40 


60 80 100 120 140 




20 40 


60 80 100 120 140 






m mo>1 [M B ] 








m mo>1 [M G ] 



Figure 8. Same as Figure[7] but for the M cl 2 simulation with M c \ i = 20Mq. 



criteria of our observed stars and clusters. For each tar- 
get star, we identify all simulated clusters having m max 
within the uncertainty of the star's empirically-derived 
mass (Table [lj . We also eliminate all simulated clus- 
ters containing more than one O star, to obtain our fi- 
nal sub-sample of simulated clusters. For 77i maXj 2, we 
use the mass of the most massive observed, compan- 
ion main-sequence star. In the case of our apparently 
isolated massive stars, we set m maXi 2 to an upper-limit 
value of 1.5Mq, our F814W completeness limit. In the 
following analysis, we denote the observed m mMi 2/™m» 
as [m maXi 2/m m ax]obs and that from simulated clusters 
as [m maXi 2/m ma x]sim- Table [7] lists the fraction of clus- 
ters with [m m ax,2/m-max] S im < [™max,2/™max]obs for each 

star. Column 1 lists the star ID; column 2 lists the 
mass of the OB star from Table [TJ column 3 lists 
[fftmax^/wimaxjobs; columns 4, 5, and 6 list the fraction 
of clusters with [r7i m ax,2/"i max ] sim < [m maXj2 /m max ) ohs 
for the simulations with N*j = 40, having slopes 
of -1.8, -2, and -2.3, respectively; and columns 7, 8, 
and 9 list these fractions for the M c \ simulations with 



M c \.\o = 20M Q , as shown. Columns 5 and 8 correspond 
to the simulated population in Figures [3> and [8^, respec- 
tively. Columns 10 and 11 are the same as Columns 5 
and 8, respectively, except that they correspond to the 
simulated population in Figures 0a and showing all 
clusters with > 1 O star. 

Two trends emerge from these data: (1) steepening 
the power law slope increases the fraction of simulated 
clusters with [m m ax, 2 /ramaxkm < [ra max , 2 /m max ] b s and 
(2) the simulations typically have a slightly lower 
fraction of simulated clusters having [m max .2 / m max ] s i m < 
[TO maXj2 /?7i ma x] bs than the M c \ simulations. The first ef- 
fect is caused by steeper slopes creating a greater fraction 
of small clusters, where stochastic effects can result in 
massive stars forming with only a few low mass compan- 
ions. The second effect is caused by the fact that some of 
the clusters fall into the low m max> 2/m max regime in the 
M c \ simulation where N* < 40, which is not allowed in 
the iV* simulation, due to our lower limit of iV*,i = 40. 
Looking specifically at the j3 = —2 simulations, Table [7] 
shows that the frequency of single O star simulated clus- 
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Table 7 

Fraction of clusters with [m maXj2 /m max ] aim < [m max , 2 /m max ] obs 



Star m m ax (M ) (^^-) obs AT" 1 ' 8 N~ 2 N~ 2 - 3 MJ 1 " 8 A/" 2 A/" 2 ' 3 7V.T 2 , > 1 Af" 2 , N Q > 1 



AzV 58 22 0.07 0.01 0.02 0.03 0.04 0.06 0.09 0.02 0.05 

AzV 67 37 0.07 0.07 0.10 0.11 0.08 0.10 0.11 0.07 0.06 

AzV 106 18 0.17 0.11 0.14 0.18 0.16 0.21 0.27 0.14 0.20 

AzV 186 33 0.05 0.01 0.02 0.02 0.03 0.04 0.05 0.01 0.03 

AzV 226 35 0.04 0.01 0.02 0.02 0.03 0.04 0.05 0.01 0.03 

AzV 302 27 0.15 0.20 0.24 0.29 0.19 0.24 0.30 0.19 0.19 



ters having [m max , 2 /m max ] s i m < [m max , 2 /m max ] b s for all 
observed groups is below the 25th percentile in both sim- 
ulations, confirming the impression from Figures [7]!) and 
[8b- For all O star clusters from Figures [TJa and [8b, our 
observed groups are all below the 20th percentile. Some 
of our isolated star observations are found in the lowest 
5th percentile in both simulations, but still within the 
parameter space covered by the simulations. 

That our observations are not well-distributed among 
the cluster population in Figures [7b and [8b is primar- 
ily due to our sample selection. Our selection process 
included a visual inspection of our targets using ground- 
based imaging to ensure they appeared isolated, thereby 
ensuring that our objects have extremely low values of 
w maX! 2/w max . These targets were drawn to qualify for 
both the field O star and field OB star samples from 
Oey et al. (2004), defined to have no other stars having 
m > 20Af© and m > IOMq within a clustering length, 
respectively. Thus by definition, rn max .2 < IOMq for our 
sample. On the other hand, 27 of the 91 stars (30%) 
in the field O star sample of Oey et al. are not mem- 
bers of the field OB sample, implying that for these stars 
w maXj 2 > 10Af Q . The remaining 70% of isolated O stars 
having companions with masses below IOMq can be com- 
pared to a simulated fraction of 57% in our M~ x 2 simula- 
tion and 64% in our iV" 2 simulation. However, the SMC 
field O star sample is contaminated by runaway stars 
which will inflate the observed fraction of O stars with 
m maX! 2 < 10M Q . 

However, we also note that the distribution of 
,2/wi max in the simulations does depends somewhat 
on the cluster population parameters and populating al- 
gorithm. For example, if we do not discard and repopu- 
late clusters when the total cluster mass exceeds 105% of 
the target mass (see £|3.2[) . then our test simulations show 
that the increased production of massive stars can affect 
the percentiles by up to a factor of two, so that our ob- 
jects fall in the lowest 50th percentile, for the default M c \ 
simulation. Thus in this case, our observed clusters are 
nearer to the median and more well-distributed overall. 

5. THE RELATION BETWEEN m max AND M cl 

5.1. Observations and Simulations 

In Figures [9] and [TOJ we compare the relationship be- 
tween the cluster mass M c \ and maximum stellar mass 
m m „. As in Figure [3 panel (a) includes simulated clus- 
ters with > 1 O star, while panel (b) includes only simu- 
lated clusters with exactly one O star. The color coding 
scheme is also the same as in Figure [7J The solid lines 
show contours for the 10th, 25th, 50th, 75th, and 90th 
percentiles of M c \ as a function of rn max in the simula- 



tion, while the dashed line represents the mean. These 
percentiles are computed from a much larger set of 10 6 
(M c i simulations) or 10 7 (AT* simulations) modeled clus- 
ters, to reduce stochastic scatter. For our observations 
(black squares), we calculate M c \ as described in i i3.ll 
(Table QJ . The diamonds show observed Galactic clus- 
ters whose M c \ and m max are tabulated by Weidner et 
al. (2010a). 

Figures IHh and fTUb show that the majority of the 
Galactic cluster sample lies above the 90th percentile in 
M c \ for a given m max in both the iV* and M c \ simula- 
tions. In contrast, our observed objects all occur below 
the 50th percentile in both simulations and are more rep- 
resentative of the single O star cluster sample in Figures 
M> and [TO]!). Indeed, Figures M> and [T0b show that the 
majority of the Galactic cluster sample lies outside the 
single O star cluster parameter space. This indicates that 
the Galactic cluster sample is comprised of clusters with 
a well-populated IMF, probably due to selection effects, 
since more fully populated clusters preferentially tend to 
be observed (Maschberger & Clarke 2008). 

In Figure [Til we plot the mean m max value as a func- 
tion of M c i (dashed line) for our M c ^ 2 simulation. The 
solid lines now show contours for the 10th, 25th, 50th, 
75th, and 90th percentiles of TO max as a function of M c \ 
in the simulation, again calculated from the larger sets of 
simulated clusters. The dashed line represents the mean. 
We note that in FigurefTTb. the percentiles arc calculated 
from all clusters having a given M c \, while the simulated 
clusters plotted here are only those which contain at least 
one O star. In FigurefTTb. the percentiles exclude all clus- 
ters with multiple O stars, and the plotted clusters are 
those with a single O star. 

Here, we see that nearly all of our observed clusters are 
above the 90th percentile of m max as a function of M c \. 
In FigurefTTb. all the Galactic clusters are below the 50th 
percentile, with the majority below the 10th percentile. 
FigurefTTb again demonstrates that the Galactic O star 
cluster sample largely falls outside the parameter space 
of the single O star clusters. Comparing the percentiles 
plotted in Figure fTTb with those in Figure fTTb shows 
little difference for our observed M c \ values; however, as 
M c \ increases, the transition from single O star clusters 
to clusters with > 1 O star is revealed in the turnover 
of the percentiles in Figure fTTb . This confirms that the 
Galactic cluster sample is comprised mainly of clusters 
with multiple O stars. 

5.2. Does M c \ determine m max ? 

We now compare our observations with numerical sim- 
ulations that are limited by a relation between the maxi- 
mum stellar mass within a cluster and total cluster mass 




Figure 9. log(m max ) vs. log(M c i) from the N* simulations for (a) all clusters having at least one O star, and (6) all clusters having only a 
single O star. Colors are as in Figure [7] with diamonds representing Galactic clusters tabulated by Weidner et al. (2010a). The solid lines 
represent the 10th, 25th, 50th, 75th, and 90th percentile of M c \ as a function of m m ax from the simulations. The dashed line represents 
the average value from the simulations. 
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(w max -M c i). As mentioned in fJTJ various forms of this 
relation have been proposed, based on both theory (Bon- 
nell et al. 2004) and observations (WK06), that invoke 
a physical relation between m max and M c \. This is dif- 
ferent from the purely statistical relation between the 
average rn max and M c \ (Oey & Clarke 2005). In the lat- 
ter case it is simply improbable to form a massive star in 
a small cluster, whereas in the integrated galaxial initial 
mass function (IGIMF) proposed by Weidner & Kroupa 
(2005), the m max -M c i relation is modeled deterministi- 
cally, such that m max never exceeds the value derived 
from this ?n max -M c i relation. For reference, see Figure 1 
of Weidner & Kroupa (2005), which plots various m max - 
M c \ relations from the literature. 

If we adopt the mean in Figure lllb as a simple mmax- 
Mci relation, then all simulated clusters above the dashed 
line are in violation of such a relation. We note that the 
mean corresponds to somewhat lower-mass clusters for 
a given m max than the WK06 m max -M c i relation. Even 
so, our observed minimal O star groups do not fit within 
the framework of a steepened IGIMF as presented by 
Weidner & Kroupa (2005), although we note that statis- 
tical variation of the m max -M c i relation is not included 
in their work. 

Fig Lire 1 121 shows ?7X maXj 2/m max a function of TTlmax 

for clusters that correspond to our imposed m max -M c i re- 
lation, which are those below the mean m max plotted in 
Figure [TTb . The color coding and panel samples are the 
same as in Figure [8] Figure fT2l a, illustrates that all sim- 
ulated O star clusters now exhibit m maXi 2/fi max > 0.2. 
With this imposed ?7i max -M c i relation, our observations 
appear to fall completely outside the parameter space 
of the simulations in both Figure [T2"b and [TTb . given 
the number of clusters corresponding to the SMC cluster 
population. This form of the deterministic m max -M c i re- 
lation is therefore poorly supported by our observations 
of minimal O star groups. 

One possible interpretation of these results is that our 
observed objects are remnants of clusters affected by "in- 
fant weight loss" (e.g., Bastian & Goodwin, 2006). How- 
ever, even if the objects have been reduced by a factor of 
a few in M c \ or N*, they would still be discrepant from 
the Galactic cluster sample of Weidner et al. (2010a). 
Recent simulations have shown that cluster mass segre- 
gation can occur on timescales of ~1 Myr for the most 
massive cluster stars (Allison et al. 2009). Massive stars 
segregated to the cluster core are unlikely to be evapo- 
rated by "infant weight loss", and so clusters are likely 
to retain the two most massive stars. Therefore, the ob- 
served m maXi 2/m max values are unlikely to be affected by 
dynamical evaporation. Moreover, "infant weight loss" 
is associated with the rapid removal of gas, which, how- 
ever, is still present in the majority of our objects ( §2.3)1 , 
Furthermore, iV-body simulations of sparse, young clus- 
ters by Weidner et al. (2010b) show that < 15% of the 
cluster's mass is removed within 5 Myr, suggesting that 
significant mass loss is relatively unimportant in such ob- 
jects. 

Our results are similar to those found by Maschberger 
& Clarke (2008), who examined a sample of clusters from 
the literature using studies that focused on low clus- 
tering around high-mass stars. The observations of iso- 
lated Herbig Ae/Be stars by Testi et al. (1997, 1999), as 
well as our observations of isolated OB stars, show that 



massive stars may form in even the most sparse environ- 
ments. As we showed above, these observations are not 
consistent with a strictly-defined, deterministic mmax- 
Mci relation. At the very least, the minimal OB groups, 
along with clusters compiled by Weidner et al. (2010a), 
imply huge deviation from a direct m max -M c i relation 
(Figure ITT]) . While a steepened IGIMF will occur even 
without a m max -M c i relation, this result has problematic 
implications for the magnitude of proposed steepening 
for the IGIMF in aggregate galactic stellar populations 
(Weidner & Kroupa 2005). In addition, the competitive 
accretion theory of star formation is also linked to an 
m max -M c i relation (Bonnell et al. 2004), although we 
note that Maschberger et al. (2010) find that compet- 
itive accretion simulations are nevertheless able to pro- 
duce massive field stars of at least 9M Q . On the other 
hand, simulations of star formation based on core accre- 
tion (e.g., Krumholz et al. 2010) show that under specific 
cloud conditions, radiative heating can prevent fragmen- 
tation, perhaps more directly forming minimal O star 
groups similar to those observed in this paper. 

6. CONCLUSIONS 

We carried out a SNAP program with HST's Advanced 
Camera for Surveys that yielded high resolution obser- 
vations of eight field OB stars in the SMC. These stars 
range in spectral type from 07 to B0. 5, and in mass from 
18M to 37M . Radial velocities for two stars (AzV 223 
and smcl6) show them to be runaways, and we may ex- 
pect one or two more to be transverse runaways. There 
is no evidence of clustering down to a 1.5M Q detection 
limit in three of the six non-runaway cases. The non- 
runaway, isolated stars (AzV 58, AzV 186, and AzV 226) 
remain candidates for isolated OB star formation. Two 
of these isolated OB stars reside within HII regions, indi- 
cating that these stars may still be located in the region 
of their formation, and strengthening the possibility that 
these O stars have formed alone. For the other three non- 
runaway OB stars (AzV 67, AzV 106, and AzV 302), we 
detect an associated population of stars using a stellar 
density analysis and a separate friends-of-friends algo- 
rithm. After accounting for field contamination, we find 
eight to ten stars associated with each OB star, ranging 
in mass from our lower detection limit of 1M Q to 4M Q . 

The three observations that do show evidence of clus- 
tering exhibit a flat IMF with slope of T = 0.1 ± 1.0 to 
— 0.2 ± 0.9 when combining their populations, although 
due to the small sample size, the IMF may be consistent 
with a Salpeter IMF, which has T = -1.35. The flat IMF 
is due to a lack of low-mass companions that ordinarily 
are expected to form along with these O stars. 

Assuming that each of our non-runaway stars is still in 
the location where it formed, we infer cluster membership 
of N,f = 19—171 based upon the companion population > 
1.5-Mq, integrated over the full stellar mass range for our 
Kroupa IMF of 0.08M Q - 150Af Q . Given their inferred 
iV», we calculate that a small fraction, only 0.01-0.2, of 
clusters will form a star with mass m > m max observed. 

We have conducted Monte Carlo simulations to explore 
where our observations fall within the cluster parame- 
ter space of typical star formation, assuming a Kroupa 
IMF and either a clustering law based on iV* or cluster 
mass function based on M c \. The power law slopes and 
lower limits of M c \ and N* for these simulations were 
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Figure 12. Same as Figure [8] but applying a fixed m max - M c \ relation. 



constrained using observations of massive star clustering 
in the SMC by Oey et al. (2004) and the fraction of iso- 
lated Galactic O stars by de Wit et al. (2005). Together, 
these observational constraints resulted in a -2 power law 
slope with lower limits of N*^ — 40 or M c \ t \ — 2OM 
as the best fit models, which we adopted as the default 
simulations for this study. These default and M c \ 
models match equally well with the observations, thus 
neither one is established as a more suitable metric for 
modeling cluster distribution. 

We find that observed mass ratios of the two highest- 
mass stars are below the 25th percentile of single O star 
clusters generated by the default simulations. This re- 
sult is due to choosing targets which appeared isolated in 
ground-based imaging. Our observations also lie below 
the 50th percentile when comparing total cluster mass 
(either M c \ or N*) as a function of m max , whereas a sam- 
ple of Galactic clusters from Weidner et al. (2010a) are 
nearly all above the 90th percentile. These numbers sug- 
gest that our observations are more typical examples of 
O star clusters than the Galactic cluster sample, which 
contains clusters with well-populated IMFs. 

We show evidence that our observed minimal O star 
groups are inconsistent with a deterministic m max -M c i 
relation. By extension, our observations are also in- 
consistent with the TOmax-Mci relation proposed for the 
IGIMF effect (Weidner & Kroupa 2005) and the relation 
that Mci oc rn^x (Bonnell et al. 2004), predicted by the 
competitive accretion model of massive star formation. 
We argue that in most cases, an observed m max -M c i re- 
lation is simply a product of the stochastic, probabilistic 
nature of a universal IMF, rather than an IMF with a 
variable upper limit corresponding to cluster mass. We 
conclude that our observations of minimal O star groups 
are consistent with a universal stellar mass function, in- 
cluding a constant stellar upper-mass limit, without the 
need to invoke a m max -M c i relation. 
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